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^ ' Abstract 

• We study the stability of a four parameter family of spatially periodic traveling wave solutions 

, of the generalized Benjamin-Bona-Mahony equation to two classes of perturbations: periodic per- 

turbations with the same periodic structure as the underlying wave, and long-wavelength localized 
perturbations. In particular, we derive necessary conditions for spectral instability to perturbations 
to both classes of perturbations by deriving appropriate asymptotic expansions of the periodic Evans 
function, and we outline a nonlinear stability theory to periodic perturbations based on variational 
, methods which effectively extends our periodic spectral stability results. 

> 
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1 Introduction 

In this paper, we consider the generahzed Benjamin-Bona-Mahony (gBBM) equation 

ut-Uxxt + Ux + if{u))^ = 0, (1) 



< 



> 

(N 
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where cc G M, t £ M, and /(•) G C'^(R) is a prescribed nonlinearity. In particular, we will be most 
interested in the case of a power-law nonlinearity f{u) = u^"*"^ as this is when our results are most 
explicit. Notice that when f{u) = u^, ^ is precisely the Benjamin-Bona-Mahony equation (BBM), or 
^\ [ the regularized long-wave equation, which arises as an alternative model to the well known Korteweg-de 

■ Vries equation (KdV) 

Ut - Uxxx + UUx = 



o 



. as a description of gravity water waves in the long- wave regime (see IHHS]). In applications, various other 

5— ( I nonlinearities arise which facilitates our consideration of the generalized BBM equation. For suitable 

nonlinearities, equation ([T]) admits traveling wave solutions of the form u{x,t) = u{x — ct) with wave 
speed c > 1 which are are either periodic or asymptotically constant. Solutions which are asymptotically 
constant are known as the solitary waves, and they correspond to either homoclinic or heteroclinic orbits 
of the traveling wave ODE 

cUxxt-{c-l)u, + {f{u))^^Q (2) 

obtained from substituting the traveling wave ansatz into ([T]) . The stability of such solutions to localized, 
i.e. L^(M), perturbations is well known [51[T51[2n]: the solitary waves form a one parameter family of 
traveling wave solutions of the gBBM which can be indexed by the wave speed c > 1. A given solitary 
wave solution Ucq {x — cot) of the is nonlincarly (orbitally) stable if the so called momentum functional 

N{c) := i /" {u,{x)^ + u'^ixf) dx 
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is an increasing function at cq, i.e. if -^Af{c)\c=co > 0, and is exponentially unstable if J^A/'(c)|c=co 
is negative. In the special case of a power-nonlinearity f{u) = u^'^^, it follows that such waves are 
nonlinearly stable if 1 < p < 4, while for p > A there exists a critical wavespeed c{p) such that the waves 
are nonlinearly unstable for 1 < c < c{p) and stable for c > c(p). 

In [l^, Pego and Weinstein found the mechanism for the instability of the solitary waves to be as 
follows: linearizing the traveling gBBM equation 



about a given solitary wave solution Uc of ([2|) and taking the Laplace transform in time yields a spectral 
problem of the form dxC[uc]v = fiv considered on the real Hilbert space L^(R). The authors then make 
a detailed analysis of the Evans function -D(/i), which plays the role of a transmission coefficient familiar 
from quantum scattering theory: in particular, D{fi) measures intersections of the unstable manifold as 
X —oo and the stable manifold as a; +oo of the traveling wave ODE. As a result, if Re(/i) > and 
D{^) — 0, it follows that the spectral problem dxC[uc]v = fiv has a non-trivial L^(R) solution, and hence 
^ belongs to the point spectrum of the linearized operatofl Using this machinery, Pego and Weinstein 
were able to prove that sign(_D(/^)) = +1 as /i ^ +oo and 



for ^ 1 and some constant C > 0. Thus, if Af{c) is decreasing at c, then by continuity there must 
exist a real /i* > such that D{fi*) = 0, which proves exponential instability of the underlying traveling 
wave. 

The focus of the present work concerns the stability properties of spatially periodic traveling wave 
solutions of ^ and in contrast to its solitary wave counterpart, relatively little is known in this context. 
Most results in the periodic case falls into one of the following two categories: spectral stability with 
respect to localized or bounded perturbations El [HI HI] , and nonlinear (orbital) stabihty with respect 
to periodic perturbations [TJ [51 [TTl [T7]. The spectral stability results rely on a detailed analysis of the 
spectrum of the linearized operator in a given Hilbert space representing the class of admissible perturba- 
tions: in our case, we will consider both corresponding to locahzed perturbations, and Lpgj.{[0, T]), 
corresponding to T-periodic perturbations where T is the period of the underlying wave. Notice that 
unstable spectrum in Lpj,j.([0, T]) represents a high-frequency instability and hence manifests itself in the 
short time instability (local well-posedness) of the underlying solution, while unstable spectrum near the 
origin corresponds to instability to long wavelength (low-frequency) perturbations, or slow modulations, 
and hence manifests itself in the long time instability (global well-posedness) of the underlying solution. 
Notice that on a mathematical level, the origin in the spectral plane is distinguished by the fact that 
the traveling wave ordinary differential equation ^ is completely integrable. Thus, the tangent space to 
the manifold of periodic traveling wave solutions can be explicitly computed, and the null space of the 
linearized operator can be built up out of a basis for this tangent space. 

Our analysis of the linearized spectral problem on Lpgj.([0, T]) parallels that of the solitary wave theory 
described above, in that we compare the large real /i behavior of the corresponding periodic Evans function 
to the local behavior near the origin, thus deriving a sufficient condition for instability. The stability 
analysis to arbitrary localized perturbations is more delicate and follows the general modulational theory 
techniques of Bronski and Johnson [8j: unlike the solitary wave case, the spectrum of the linearized 
operator about a periodic wave has purely continuous spectrum and hence any spectral instability 
must come from the essential spectrum. As a result, there are relatively few results in this case. In the well 
known work of Gardner 12J, it is shown that periodic traveling wave solutions of ([T]) of sufRciently long 
wavelength are exponentially unstable whenever the limiting homoclinic orbit (solitary wave) is unstable. 

^By a standard argument, the essential spectrum can bo shown to he on the imaginary axis, and hence any spectral 
instability must come from the discrete spectrum. 



Ut - Uxxt + CU. 



(c - l)Ux + f{u)x = 



(3) 
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The mechanism behind this instabihty is the existence of a "loop" of spectrum in the neighborhood of any 
unstable eigenvalue of the limiting solitary wave. More recently, Haragu§ carried out a detailed spectral 
stability analysis in the case of power- nonlinearity f{u) = uP~^^/{p+ 1) for waves sufficiently close to the 
constant state u — ((p + l)(c— l))""^^^ determining when the spectrum of the linearized operator is 
confined to the imaginary axis. We, on the other hand, consider arbitrary periodic traveling waves with 
essentially arbitrary nonlinearity: as a result of this level of generality we are unable to make conclusive 
spectral stability statements but instead can only determine L^-spectral stability near the origin: this 
is determined by an asymptotic analysis of the periodic Evans function near the origin resulting in a 
modulational instability index which determines the local normal form of the spectrum. However, we 
find that this analysis yields quite a bit of information about the spectrum of the underlying wave. 

The orbital stability results rely on the now familiar energy functional techniques or Grillakis, Shatah, 
and Strauss [13] used in the study of solitary type solutions. This program has recently been carried out 
by Johnson [T7] in the case of periodic traveling waves of the generalized Korteweg-de Vries equation. 
The corresponding theory for the gBBM equation is nearly identical to the analysis in [T7], and hence 
will only be outlined in this work. 

Finally, we study the behavior of the above indices in a long- wavelength limit in the case of a power- 
nonlinearity f(u) — u^'^^ /{p + 1). For such non-linearities, we gain an additional scaling in the wave 
speed which allows explicit calculations of the leading order terms in this stability index. A naive guess 
would be that the value of the orientation index would converge to the solitary wave stability index. This 
is not the case, however, since the convergence of long- wavelength periodic waves to solitary waves is non- 
uniform, implying such a limit is highly singular. What is true is that the sign of the finite-wavelength 
instability index converges to the sign of the solitary wave stability index. Since this index was derived 
from an orientation index calculation, only its sign matters and thus we show that this stability index 
is in some sense the correct generalization of the stability index used in the study of solitary waves of 
([1]). Moreover, the sign of the modulational stability index converges to the same sign as the solitary 
wave index, implying that periodic waves of ([T]) in a neighborhood of a solitary wave are modulationally 
unstable if and only if the nearby solitary wave unstable. Notice this does not follow directly from 
the work of Gardner described above: there it is proved that there exists a "loop" of spectrum in the 
neighborhood of the origin which is the it continuous image of a circle, but to our knowledge it has never 
been proved that this map is injective. 

The outline for this paper is as follows. In section 2 we review the basic properties of the periodic 
traveling wave solutions of ([l]), and in section 3 we review the basic properties of the periodic Evans 
function utilized throughout this work. In section 4, we begin our analysis by considering stability of the 
T-periodic traveling wave to T-periodic perturbations. In particular, we first determine the orientation 
index, which provides sufficient information for spectral instability to such perturbations, and then discuss 
how this index plays into the nonlinear stability theory. In section 5, we conduct our modulational 
instability analysis. In section 6 we analyze the results of sections 4 and 5 in a solitary wave limit, thus 
extending the well known results of Garnder. Finally, we close in section 7 with a brief discussion and 
closing remarks. 

2 Properties of the Periodic Traveling Waves 

In this section, we describe the basic properties of the periodic traveling waves of the gBBM equation 
([T]). For each c > 1, a traveling wave is a solution of the traveling wave ODE ([2]), i.e. they are stationary 
solutions of ([I]) in a moving coordinate frame defined by a; — ct. Clearly ^ defines a Hamiltonian ODE 



Stability of Periodic gBBM Waves 



4 



and can be reduced to quadrature: in particular, the traveling waves satisfy the relations 

CMxa; " (c - l)u + f (u) = fl 

lul-(^^yi'^ + Fiu) = au + E (4) 

where a and E are real constants of integration and F' ^ f with F{Q) = 0. Notice that the solitary waves 
correspond to a = and E fixed by the asymptotic values of the solution. However, the parameters a 
and E are free parameters: we must only require that the effective potential 

c — 1 

V{u; a, c) :— F{u) — — au 

have a non-degenerate local minimum (see Figure [T]). Note that this places restrictions on the allowable 
parameter regime for our problem: we will always assume we are in the interior of this region, and that 
the roots u± of the equation V{x;a,c) — E are simple and such that V{x;a,c) < E for x € 
In particular, this guarantees the classical "turning points" u± are functions of the parameters a, E, 
and c. Thus, the periodic solutions of ^ form a four parameter family of solutions u(x + xq] a, E, c) 
while the solitary waves form a codimension two subset. Notice however that the translation invariance, 
corresponding to the parameter xq is not essential to our theory and can be modded out. Hence, we 
consider the periodic traveling wave solutions of ([T]) as a three parameter family of the form u{x; a, E, c). 
The partial differential equation ([1]) has, in general, the three conserved quantities 

M = / {u — Uxx) dx 
Jo 

P = II {u' + ul)dx (5) 



which correspond to the mass, momentum, and Hamiltonian (energy) of the solution, respectively. These 
three quantities are considered as functions of the traveling wave parameters a, E, and c and their 
gradients with respect to these parameters will play an important role in the foregoing analysis. It is 
important to notice that when restricted to the four-parameter family of periodic traveling wave solutions 
of (H]), the mass can be represented as M = u dx. Since all our results concern this four-parameter 
family, we will always work with this simplified expression for the mass functional. 

As is standard, one can use equation ^ to express the period of the periodic wave u as 



T{a,E,c) ^2^c I 

J U- 



du 



v/2 {E~V{u;a,c)) 



The above interval can be regularized at the square root branch points u± by the a standard procedure 
(see PI for example) and hence represents a function of (a, E, c). Similarly, the mass and momentum 
can be expressed as 



M{a,E,c) = / 

J U 



"+ u du 



v/2 {E-V{u;a,c)) 



and can be regularized as above. In particular, it follows that one can differentiate these functionals 
restricted to the periodic wave u{x;a,E,c) with respect to the parameters {a,E,c). The gradients of 
these quantities will play an important role in the subsequent theory. 
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Figure 1: (Left) A plot of the effective potential energy V{x;Q.l,2) for the modified BBM equation 
(/(w) — u^), as well as three energy levels Ei = —0.32, E2 = —0.05, and E3 = 0.1. (Right) Plots in phase 
space of the solutions u{x;0.1, Ej,2) corresponding to the three energy levels on the left. Notice that 
those solutions corresponding to energy levels Ei and E2 are bounded by a homoclinic orbit in phase 
space (given by the thin dashed line). Moreover, notice that E2 corresponds to two distinct periodic 
traveling wave solutions: however, these can be clearly distinguished by their initial values which we have 
chosen to mod out in our theory. 



It is useful to notice the following connection to the classical mechanics corresponding to the traveling 
wave ODE. The classical action in the sense of action angle variables is given by 

2 



K{a,E,c) = f u^du ^ 2J - y/{E - V{u;a,c))du. (6) 



While K is not itself conserved, it does provide a useful generating function for the conserved quantities 
of ID). Indeed, it is clear the classical action satisfies the relation 

Va,E,cKia, E, c) = ( -M (a, E, c), -T(a, E,c),- (P(a, E, c) - K{a, E, c)) 
\ c c c 

where V a,E,c {da, ds, dc) , which immediately establishes several useful identities between the gradients 
of T, M, and P. For example, it follows that Ta — Me and Tc = Pe- 

Finally, we make a few notes on notation. Throughout the forthcoming analysis, various Jacobians of 
maps from the traveling wave parameters to the period and conserved quantities of the gBBM flow will 
become important. We adopt the following Possion bracket style notation 

r ,x 9{g,h) 
^^'"^--^ 

for the Jacobian determinants with the analogous notation for larger determinants 

r , .1 d{g,h,j) 



Notice that the nonvanishing of such quantities encodes geometric information about the underlying 
manifold of periodic traveling wave solutions of ((T]): more will be said on this in the coming sections. 
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3 The Periodic Evans Function 

We now begin our stability analysis of a T = T{a, E, c)-pcriodic traveling wave u{x; a, E, c) by considering 
a solution to the partial differential equation i.e. a stationary solution of ([3]), of the form 

tp{x, t) = u{x; a, E, c) + ev{x, t) + O(e^) 

where jej ^ 1 is considered as a small perturbation parameter. Substituting this into ([1]) and collecting 
terms at 0{e) yields the hnearized equation JC[u]v — Vvt where J = dx, C[u] — —cd^ + (c — 1) — f'{u), 
and D = 1 — 9^. Since this linearized equation is autonomous in time, we may seek separated solutions 
of the form v{x,t) = e'^*v{x), which yields the spectral problem 

J£v = fiVv (7) 

Throughout this paper, we consider the above operators as acting on _L^(R) corresponding to spatially 
localized perturbations, or on Lpj,^([0,T]) corresponding to T-periodic perturbations. In both cases, the 
operator I? is a positive operator and is hence invertible and hence ([7]) can be written as a spectral 
problem for a linear operator: 

Av = fiv, A = V^^JC. 

As £ has is a differential operator with periodic coefficients, the natural setting to study the L^(M) 
spectrum of the operator A is that of Floquet theory. We begin with the following standard definition. 

Definition 1. The monodromy operator M(fj,) is defined to be the period map 

M{^l)^^T■,^i) 

where fi) satisfies the first order system 

<i>{x;n)x = U{x,^)<i>{x;n) (8) 
subject to the initial condition $(0) = I, where I is the 3x3 identity matrix and 

/ 10 

ll{x,fi) = 1 

V + i(c-l-/'(^.)) f 

It now follows from an easy calculation that ([7]) has no point spectrum in L^(R). Indeed, suppose 
^ is a vector solution of ([8]) corresponding to a non-trivial L^(R) eigenfunction of A with eigenvalue /i. 
From the definition of the monodromy operator, we have 

*(A^T) = M(^)^*(0) 

for any N E Z. It follows that ^ can be at most bounded on R and must not decay as ^ ±00. This 
observation leads one to the following definition. 

Definition 2. We say /i € spec(y^) if there exists a non-trivial bounded function such that Aip = fJ-tp 
or, equivalently, if there exists a X such that 

det(M(^)-AI) = 0. 

Following Gardner [IB^ we define the periodic Evans function _D : C x C — > C to be 

D{fi,X) ==dct(M(/^)- AI). 

Finally, we say the periodic solution u(x; a, E, c) is spectrally stable if spec{A) does not intersect the open 
right half plane. 
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Remark 1. First, notice by the Hamiltonian nature of ([7]), spec(y^) is symmetric with respect to reflec- 
tions about the real and imaginary axis. Thus, spectral stability occurs if and only i/spec(^) C Mi. 

Secondly, since we are interested primarily with the roots of D[ix, A) for A on the unit circle, we 
will frequently work with the function D{iJ,,e^'^) for k G M, which is actually the function considered by 
Gardner. 

It follows that we can parameterize the continuous (R) spectrum of the operator A by the Floquet 
parameter k: 

spec{A) = U {^IeC■. D{n, e'"") = 0}. 

K(E [ — TT, 7r) 

In particular, the zero's of the function _D(/i, e*'') for a fixed k g R correspond to the Lp^^{[0,T])- 
eigenvalues of the operator resulting from the map dx i-^ dx+iK applied to ([7]), and hence the continuous 
spectrum of A can be parameterized by a one-parameter family of eigenvalue problems. By a standard 
result of Gardner, if D(/Lto, e*"") = 0, then the multiplicity of fiQ as a periodic eigenvalue of the corre- 
sponding linear operator is precisely the multiplicity of fiQ as a root of the Evans function. As we will 
see below, the integrable structure of ^ implies that the function D{fi, 1) has a zero of multiplicity 
(generically) three at /i = 0. For small k then, there will be in general three branches of roots of 

D{fij{K), e*") which bifurcate from the origin. Assuming these branches are analytic^ in k, it follows that 
a necessary condition for spectral stability is thus 

^M,(«)|_„GMz. (9) 

This naturally leads to the use of perturbation methods in the study of the spectrum of A near the origin, 
i.e. modulational instability analysis of the underlying traveling wave. As we will see, the first order terms 
of a Taylor series expansion of the three branches /^^(k) can be encoded as roots of a cubic polynomial, 
and hence spectral stability is determined by the sign of the associated discriminant. Moreover, it follows 
by the Hamiltonian structure of ^ that in fact a{A) C Mi if (O holds and the roots of the cubic 
polynomial are distinct. 

We conclude this section by reviewing some basic global features of the spectrum of the linearized 
operator A which are useful in a local analysis near fj, = 0. We also state some important properties of 
the Evans function A) which are vital to the foregoing analysis. 

Proposition 1. The L'^{M.)- spectrum of the operator A has the following properties: 

(i) There are no isolated points of the spectrum. In particular, the spectrum generically consists of 
piecewise smooth arcs. 

(ii) The entire imaginary axis is contained in the spectrum, i.e. iR C spec(yl). 

Moreover, the Evans function D{p,,X) satisfies the following: 

(lit) D{fi,0) = e^^/^ 

(iv) D{^i, A) ^ det( M(^) - AI) = -A^ + a{^i)\^ - a(-^i)e^^/=A + e^^'^'" with a{fi) = tr( M(^)). 

Proof. The first claim, that the spectrum is never discrete, follows from a basic lemma in the theory of 
several complex variables: namely that, if for fixed A* the function D^n, A*) has a zero of order k at /i* 
and is holomorphic in a polydisc about (/i*. A*) then there is some smaller polydisc about (/i*. A*) so that 
for every A in a disc about A* the fmiction A) (with A fixed) has k roots in the disc \fi — /^*| < S. 

^In general, for each j, the theory of branching solutions of non-linear equations guarantees the existence of a natural 
number nij such that /ij(-) is an analytic function of k^/^j . As we will see in our case, the Hamiltonian nature of the 
linearized operator A assures that mj = 1, and hence the roots are in fact analytic functions of the Floquet parameter. 
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For details see the text of Gunning[14|. It is clear from the implicit function theorem that ^ is a smooth 
function of A as long as ^ = tr(cof*( M(/x) — AI) M|J 7^ 0, where cof represents the standard cofactor 
matrix. 

Claim [in) follows from Abel's formula and the fact that ctr( H(a;, /i)) — /i, where H(a;,/z) is as in 

©• 

Next, we prove claim (iv). Since Av = iiv is invariant under the transformation x —x and ^ 1-^ — /i, 
we have M{ii) ~ M{—ij)~^. If we define a(/i) as above and b{ii) such that 

det[M(^) - AI] = -A^ + a{^l)\^ + 5(^)A + e^^^", 

it follows that 

det[M(^)-AI] = det[M-\-/i)-AI] 

= -A'''det[M^\-/i)]det[M(-/x) - A-i] 

= -A^ - e^^/^6(-/i)A2 - e''^/^a(-/i)A + e'^^/^ 

Therefore, b{fi) = — e^"^/'^a(— /i) as claimed. 

Claim (m) now follows from a symmetry argument. Since a((u) is real on the real axis, it follows by 
Schwarz reflection that for /i G we have a(/l) = a(/i). For /i e Mi then the Evans function takes the 
form 

D{fi, A) = -A^ + a(/i)A2 - e^'^^A + e'''^/'=. 
It follows that 

so that the roots of -D(/i, A) for a fixed /i e Ri are symmetric about the unit circle. Since there must be 
three such roots, one must live on the unit circle and hence fi £ spec(yl) as claimed. □ 



4 Periodic Instabilities: Spectral and Nonlinear Stability Re- 
sults 

In this section, we make a detailed analysis of the stability of a given periodic traveling wave solution 
of ll]) to perturbations with the same periodic structure: such perturbations correspond to A = 1 in the 
above theory. We begin by a spectral stability analysis, and then conclude with a brief discussion of a 
nonlinear (orbital) stability result. 



4.1 Periodic Spectral Instabilities 

Let a, E, c) be a T-periodic traveling wave solution of ([T]). Considering the spectral stability of such 
a solution of perturbations which are T-periodic is equivalent to studying the spectrum of the linear 
operator A on the real Hilbert space Lp^^.{[0, T]). We begin with the following lemma. 

Lemma 1. Let u{x; a, E, c) be the solution of the traveling wave equation ^ satisfying u(0; a, E, c) = w_ 
and Ux{0] a, E, c) = 0. A basis of solutions to the first order system 

= H(x,0)y 



cux{x;a, E,c) \ / cua{x]a, E,c) \ / cuEix;a, E,c) 

cUa;x{x;a,E,c) , 1^2(2;)= cUax{x;a,E,c) , y3(x) = cuEx{x;a,E,c) 
cuxxx{x;a,E,c) J \ cUaxx{x;a,E,c) J \ cuExx{x;a,E,c) 
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Moreover, a particular solution to the inhomogeneous problem 

= H(x,o)y + VF 

where W — {Q^Q,c'Dux) is given by 

(-cuc{x; a, E, c) 
-~cucx{x; a, E, c) 
-cUcxx{x;a,E,c) 

Proof. This is easily verified by differentiating ^ with respect x and the parameters a, E, and c. □ 
By the above lemma, three linearly independent solutions of the differential equation 

dxCv = 

are given by the functions u^^ Ua, and ue above. Notice that here we are considering the formal operators 
with out any reference to boundary conditions. In order to understand the structure of the periodic null- 
space of the operator we now use Lemma [1] to express the solution matrix in this basis at a; = T and 
a: = 0. 

Notice by hypothesis, for any a, E, c the solution u satisfies 

u{0;a,E,c) = ^ u{T]a,E,c) (10) 
Ux{0;a,E,c) = ^ Ux{T,a, E,c) (11) 

Uxx{0;a,E,c) = --V'{u^-a,c)^Uxx{0;a,E,c) (12) 
c 

and, moreover, it follows from ^ that Uxxx{0', a, E, c) = 0. Defining U(a:, 0) = [Yi (x), Y2{x), i1j(x)] where 
Yi, Y2, Y3 are vector functions corresponding to the solutions Ux, Ua, and ue, respectively, it follows by 
differentiating the above relations that 

/ c— c— 

da ^ SB 

U(0,0) = -V'{u-) 1 . (13) 

V l + (c-l-/'(^_))^ -/'(«_)) ^ 

Differentiating the relation E = V'{u^:a,c) with respect to E gives 

9u_ 



det(U(0,0)) = -cV'{u-)' 



dE 



and hence these solutions are linearly independent at x = 0, and hence for all x. Moreover, using the 
chain rule to differentiate pHll - p^ . we see the matrix U(r; 0) is given by 



U(T,G) = U(0,0)+ I V'{u_)Ta V'{u.)Te 


It follows that U(T, 0) — U(0, 0) is a rank one matrix. 

Our immediate goal is to relate this information to the structure of the periodic Evans function 
Z?(/x, e'"). As we will see, the fact that M(0) — I is rank one allows for significant simplifications in 
the perturbation calculations. Without this fact, one would have to compute variations in the vector 
solutions at = to an extra degree, which would force the use of multiple applications of variation of 
parameters along the same null direction. In particular, we have the following lemma. 
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Lemma 2. For \fj.\ <C 1, the periodic Evans function satisfies 

D{ii, 1) - -{T, M, P}a,£,c/^' + O(ImI')- 

Proof. Let Wi{x;fj,), i = 1,2,3, be three linearly independent solutions of ([H]), and let W(x,/i) be the 
solution matrix with columns Wi. Expanding the above solutions in powers of fi as 

w^{x, fi) = w°{x) + nwlix) + fi^wf{x) +0{\n\^) 

and substituting them into ((H]), the leading order equation becomes 

^wUx) ^ U{x,0)wnx). 
Using Proposition [l] we choose w^(a;) = yi(a::).The higher order terms in the above expansion yield 

^wi(x) - n{x,0)wi{x) + Vr\x), J > 1, (14) 

where V^"'^^ = ^0, 0, — c~^I?(u)^^^)i^ and {v)i denotes the first component of the vector v. Notice that 

for each of the higher order terms j > 1, we require (0) = 0. Notice this implies that W(0, /i) = U(0, 0) 
in a neighborhood of /i = 0, where U(0, 0) is defined in (jlSp . The solution of the inhomogeneous problem 
is given by the variation of parameters formula 



wlix) = W(x,0) 



/ w{s,oy^vr\s)ds 

Jo 



cux 'D{wf~^)i{u,u.j:}a^E dz - Ua 'D{wf~^)idz + Ue 'D{wl^^)iu dz \ 

CUxx Jq 'D{wj~^)l{u,Ux}a,E dz - Uax 'D{wl~^)idz + UEx Jq T>{wj~^)iu dz (15) 
CUxxx Jq 'D{'wl~^)i{u,Ux}a,E dz - Uaxx 'D{'wj~^)idz + UExx T>{w^.^^)iu dz / 

for j > 1. Notice we have used the identities c{u,Ux}e,x = —1 and c{u,Ux}x,a ~ m extensively in the 
above formula, which can be easily derived via equation Indeed, differentiating ^ with respect to 
E and subtracting ueUxx immediately yields the first identity. 

Now, notice that it would be a daunting task to use to the specified order needed. However, the 
integrable structure of ([7]) allows for an alternative, yet equivalent, expression in the case i ~ j ~ 1 which 
makes a seemingly second order calculation come in at first order. Indeed, in this case equation (fT4|) is 
equivalent to Lf)w\ ~ Ux and hence it follows from Proposition [T] that we can choose 

/ -CUc \ . ^ \ / \ 1 \ / 

w\{x)^ \ ^cu', + U_ + -F'(w_;a,c) j c< - ( — + -F'(m_; a, c) j 
\ / \ cu^ I 

where the above constants in front of Y2 and I3 are determined by the requirement Wi(0) ~ 0. Thus, one 
can determine the second order variation of w\ in [i by using \Vb\ to compute the first order variation of 
the function w\ defined above. Defining (5W(/i) :— W(x, /i)|^_jj, it follows that W(/i) can be expanded 
as 

0{i?) 0(m) 0(/i) 

lxV'{uJ)P{u-) + 0{i?) V'{uJ)Ta V'{uJ)Te 
O(m') 0{ii) 0(/i) 



where P{x) = —T^ + (x + ^^/'(a;)) Ta — \^ + \y' {xyj Te and the higher order terms are determined by 
(HH). Thus, 

^(/x, 1) = det ((5 W(/z) W(0, 0)-i) = 0{iX^) 



Stability of Periodic gBBM Waves 



11 



and in particular a straightforward calculation yields 

det {SW{fi)) = c{T,M,P}a,E,c + 0{fi^). 

The proof is complete by recalling that det( W(0, 0)) = — c. □ 

Remark 2. Notice that the formula for D^f^^{0, 1) differs from that derived by Bronski and Johnson 
in the case of the generalized KdV by a factor of one-half which comes from that fact the differences in 
the definitions of the momentum functionals in each work: in particular, the factor of ^ in ^ is not 
present in 

Assuming the Jacobian {T, M, P}a,E,c is non-zero, it follows that zero is a T-periodic eigenvalue of A 
of multiplicity three. By computing the orientation index 

sign ({T, M, P}a,E,c) sign (i?(oo, 1)) 

it is clear that we will have exponential periodic instability of the underlying wave if this index is negative. 
With this in mind, we now compute the large /i 3> 1 behavior of the periodic Evans function in the next 
lemma. 

Lemma 3. The function I?(-,l)|jg:M— s-K satisfies the following asymptotic relation: 

lim -Df/i, 1) — —oo. 

Proof. This follows from a simple calculation. By introducing variable change y = |/x|~-'^/^a; in ([7|, it 
follows that as /i ^ +oo the monodromy operator satisfies the asymptotic relation 

M(M)~exp(A(M)T/c), |Ai|»l, (16) 

where 

/ c 

A{n) = c 
\ -^1 n 

The eigenvalues Xj of A{ii) are rather complicated, but for /i large and real they satisfy Re(Aj) ~ a-jfJ-, 
where 

fli = 1 

1 1 21/3 

02 = as 



3 322/321/3 621/3' 
It follows that 

and thus, since 02 = 03 < < ai < c, the leading order term of 

D{^l, 1) = -1 + a{n) - a(-Ai)e^^/^ + e'^^/^ 

as /i — * +00 comes from — a(— /i)e^"^/'^, which completes the proof. □ 

We are now able to state our first main theorem for this section. 

Theorem 1. let u{x] oq, Eq, cq) be aperiodic solution of If {T, M, P}a^E,c is negative at (oq, Eq, cq) , 
then the number of roots of D{fiA) ft-^- the number of periodic eigenvalues of A) on the positive 
real axis is odd. In particular, if the period is an increasing function of energy at (gq, Eq,co), i.e. 
TE{ao, Eq, Co) > 0, then the periodic traveling wave u{x; oq, Eq, cq) is spectrally unstable to T{ao, Eq, cq)- 
periodic perturbations if and only if {T, M, P}a,E,c is negative at {clq, Eq,cq) . 
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Proof. By our work in the proof of Lcmma[21 we know that D{fi, 1) ~ Od/ip). Thus, if D^f^f^{0, 1) < 
for small positive n, the number 1) is negative for small positive fi. Since D{^, 1) is positive for real 
fi sufficiently large, we know that D{zLiJ,* , 1) = for some fi* G M*. Moreover, since Cu^ = it follows 
that C has either one or two negative T-periodic eigenvalues. Since the number of negative eigenvalues of 
C provide an upper bound for the number of unstable T-periodic eigenvalues of A with positive real part 
(see Theorem 3.1 of [H]), and since Te > implies £ has precisely one negative T-periodic eigenvalue 
(see Lemma 4.1 of [T7]) the proof is now complete. □ 

We now wish to give insight into the meaning of {T, M, P}a.E.c = at the level of the linearized 
operator A. To this end, we consider the linearized operator A as acting on Lpg^(0, T), the space of 
T-periodic functions on R. To begin, we make the assumption that {T, M}a,E and {T, P}a,E do not 
simultaneously vanish. This assumption will be shown equivalent with the periodic null-space reflecting 
the Jordan structure of the monodromy at the origin. Throughout this brief discussion, we assume that 
{T, M}a.E 0: trivial modifications are needed if {T, M}a.E vanishes but {T, P}a,E does not. First, 
define the functions 

(^0 = {T, u}a,E, "00 = 1, 

^ {T,M}a,E Uj;, V'l = / T>(t)2{s)ds, 

Jq 

4^2 = {u,T,M}a,E,c ^2 = -{T,M}E.,c+{T,M}a,EVu, 

Clearly, each of these functions belong to Tp(,j.([0, T]) and 

A(j)o = AHo = 

A<l)i = ^t^i = '02 

A4>2 = -4)1 AH2 = 0. 

In particular, we have used the fact that I?~^(l) = 1 on Lpg^([0,T]). Thus, it follows that the periodic 
null space of ^ is generated spanned by the functions 4>o and 4>i- Moreover, since 

{Mo)={T,M}a.E 
{40,'DU) ^{T,P}a^E 

the assumption that {T, M}a.E and {T, P}a,E do not simultaneously vanish implies that iVpcr {A"^) — 
Nper (A) = span((/)2), thus reflecting the Jordan normal form of the period map at = 0. 

Finally, we study the structure of the generalized periodic null space, and seek conditions for which 
there is no non-trivial Jordan chain of length two. By the Frcdholm alternative, such a chain exists if 
and only if 

{Vu,42)^{T,M,P}a.E.c^0. 

Thus, the vanishing of {T, M, P}a.E.c is equivalent with a change in the generalized periodic-null space of 
the linearized operator A. This insight has a nice relationship with formal Whitham modulation theory. 
One of the big ideas in Whitham theory is to locally parameterize the periodic traveling wave solution 
by the constants of motion for the PDE evolution. The non-vanishing of certain Jacobians is precisely 
what allows one to do this. In fact, the non-vanishing of {T, M, P}a,E,c is equivalent to demanding that, 
locally, the map (a, E, c) i— > (T, M, P) have a unique inverse: In other words, the constants of motion 
for the gBBM flow are good local coordinates for the three-dimensional manifold of periodic traveling 
wave solutions (up to translation). Similarly, non- vanishing of {T, Afja^g and {T, P}a,B is equivalent to 
demanding that the matrix 

/ Ta Ma Pa \ 

\ Te Me Pe J 
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have full rank, which is equivalent to demanding that the map (a, E) t-^ (T, M, P) (for fixed c) have a 
unique inverse, i.e. two of the conserved quantities provide a smooth parametrization of the family 
of periodic traveling waves of fixed wave-speed. 

To summarize, the vanishing of {T, M, P}a.E.c, is connected with a change in the Jordan structure of 
the hnearized operator A considered on L'^^^[[Q, T]). In particular, {T, M, P}a.E.c ensures the existence of 
a non-vanishing Jordan piece in the generalized periodic null-space of dimension exactly one. Moreover, 
and perhaps more importantly, it guarantees that infinitesimal variations in the constants arising from 
reducing the family of periodic traveling waves to quadrature are enough to generate the entire generalized 
periodic null-space of the linearized operator A: such a condition is obviously not necessary in our 
calculations, but provides significant simplifications in the theory. 

4.2 Nonlinear Periodic-Stability 

We now compliment Theorem [T] by considering in what sense the Jacobian {T, Af , P}a,E,c affects the 
nonlinear stability of a periodic traveling wave solution of ^ to T-periodic perturbations. Clearly, 
the positivity of this index is necessary for such stability but it is not clear if this is also necessary. 
Indeed, the analysis presented below allows for the possibility that a periodic traveling wave which is 
spectrally stable to T-periodic perturbations could be nonlinearly unstable to such perturbations: such a 
result would stand in stark contrast to the solitary wave theory where these two notions of stability are 
equivalent (assuming the nondegeneracy condition ^A/'(c) ^ 0). More will be said on this at the end of 
this section. 

Such analysis has recently been carried out in the context of the generalized Korteweg-de Vries equa- 
tion by Johnson |17j . There, sufficient conditions for nonlinear stability to T-periodic perturbations were 
derived in terms of Jacobians of various maps between the parameter space (a, i?, c) and the period, mass, 
and momentum. As the theory for the gBBM equation is nearly identical to this work, we only review 
the main points of the analysis here and refer the reader to [17) for details. 

To begin, we assume the nonlinearity / is such that the Cauchy problem for ^ is globally well 
posed on a real Hilbert space X of T-periodic functions defined on R, which we equip with the standard 
T^([0,T]) inner product. In particular, we require Tppj.([0,T]) to be a subspace of X. Also, we identify 
the dual space X* through the usual pairing. Now, we fix a periodic traveling wave solution of ([T]) 
uq{x; oq, -Eq, Co) and notice that we can write the linearized spectral problem ([7]) about uq in the form 

where fo(0) is an augmented energy functional defined by 

£o(0) = -^(0) + coV{(b) + aoM{<j)) 
where £, V, and M are the energy, momentum, and energy functionals defined on X defined by 

^^4')= [ {]^u{xf + F{u{x))^ dx 

M{<j)) ^ / u{x)dx 
Jo 

^(0) = ^^ {(f>{xf+M^f)dx. 

In particular, notice that each of these functionals are left invariant under spatial translations and hence 
it is appropriate to study the stability of periodic solutions of ^ up to translation. To this end, we 
introduce a semi-distance p : X x X ^ M. defined via 

p{^,lb) := inf ||0(.)-^(-+C)|U 
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and we seek conditions for which the foUowing statement is true: If Cz X is near uq as measured by 
the semidistance p, then the solution 0(a;, t) of ^ with initial data 4>o stays close to a translate of uq for 
all time. 

We begin by noticing that uq{x; oq, i?o, cq) is a critical point of the augmented energy functional Sq. 
In order to determine the nature of this critical point, it is necessary to analyze the second derivative 
of Eq. If £q{uo) £ £{X,X*) is positive definite, then nonlinear stability follows by standard arguments. 
However, one sees after an easy calculation that 

which is clearly not positive definite by the translation invariance of ([T]). Indeed, since Cuq^x = and uq 
is not monotone it follows by standard Sturm-Liouville arguments that zero is either the second or third 
eigenvalue of £ with respect to the natural ordering on M. By Lemma 4.1 of [MJl], we it follows that 
C considered on Lp^^{[0,T]) has precisely one negative eigenvalue, a simple eigenvalue at zero, and the 
rest of the spectrum is positive and bounded away from zero if Tg > 0. If Te < 0, then either the the 
null space or the number or negative eigenvalues jumps by one: a situation which can seemingly not be 
handled by the present variational analysi^. 

Thus, assuming Te > it follows that uq is a degenerate critical point of £q with one unstable 
direction and one neutral direction. In order to get rid of the unstable direction, simply notice that the 
evolution of ([1]) does not occur on the entire space X, but on the codimension two subset 

^0 ■■= {(t> e X : M{^) = M(ao, £^0, co), V{(t>) = P{ao, Eo, cq)}. 

Clearly Eq is a smooth submanifold of X containing all translates of the function uq. Defining Tq to be 
the tangent space of Sq at uq, i.e. 

To :={0eX: (uo,0> = (1,<^> -0}, 

we have by Lemma 4.3 of [MJl] that the quadratic form induced by £q{uq) is positive definite on Tq if 
Te > and 

(5^'(mo)02, 02> = -{T, M}a.E{T, M, P}a,E,c < (17) 

where (j)2 is defined as above. Since the underlying periodic wave is spectrally unstable if {T, Af , P}a,E,c < 
by TheoremlU the only interesting case in which (fT7|) holds is when {T, M}a.E and {T, M, P}a,E,c are 
positive. With these conditions in mind, it follows that the augmented energy £q is coercive on ctq near 
Uq with respect to the semi-distance p, from which nonlinear stability follows: see Proposition 4.2 and 
the corresponding proof in [17 for details. Summarizing, we have the following orbital stability result. 

Theorem 2. Let u{x; ao, Eq, cq) be a periodic solution of ^ such that the Jacobians Te, {T, M}a^E, 
and {T, M, P}a^E,c are positive at {ao, Eo,co)- Then there exists positive constants Cq and Sq such that 
if £ X satisfies p{(j)o,u) < e for some e < Eq, then the solution (j){x,t) of ((TJ with initial data 4>o 
satisfies p{(l){-,t),u) < CqE. 

By Theorem [21 it follows that {T, M, P}a.E.c > may not be sufficient for nonlinear stability of a 
periodic traveling wave of llj. In particular, notice that in the solitary wave theory it is always true that 
the operator C has only one L^(R) eigenvalue, and hence it is always true that any unstable eigenvalues of 
the linearized operator must be real. In the periodic context, however, we see this is only true if Te > 0, 
which is not true for all periodic traveling wave solutions of llJ: for example, it is clear that the cnoidal 
wave solutions of the modified BBM equation corresponding to ((U with f{u) = of sufficiently long 
wave period satisfy Te < 0. Thus, it may be possible in certain situations that the linearized operator 
A has unstable T-periodic eigenvalues which are not real. Moreover, even if Te and {T, M, P}a,E,c are 

^However, see the recent work of Bronski, Johnson, and Kapitula [BrJK]. 
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positive, and hence one has periodic spectral stabihty, it is not clear from this analysis whether one has 
nonlinear stability since the sign of the Jacobian {T, M}a^E still plays a seemingly large rol^. This 
phenomenon, which stands in contrast to the solitary wave theory, is a reflection of how the periodic 
traveling wave solutions of ([T]) have a much richer structure than the solitary waves, allowing for possibly 
more interesting dynamics. 

5 Modulational Instability Analysis 

In this section, we begin our study of the spectral stability of periodic traveling wave solutions of the gBBM 
equation ([1]) to arbitrary localized perturbations. In particular, our methods will detect instabilities of 
such solutions to long wavelength perturbations, i.e. to slow modulations of the underlying wave. Such 
stability analysis seems to be a bit more physical than the periodic stability analysis conducted in the 
previous section in the following sense: in physical applications, one should probably never expect to find 
an exact spatially periodic wave. Instead, what one often sees is a solution which on small space-time 
scales seems to exhibit periodic behavior, but over larger scales is clearly seen not to be periodic due to 
a slow variations in the amplitude, frequency, etc, i.e. one sees slow modulations physical parameters 
defining the solution. Thus, it seems natural to study the stability of such solutions by idealizing them as 
exact spatially periodic waves and then study the stability of this idealized wave to slow modulations in 
the underlying parameters. This is precisely the goal of the modulational stability analysis in this paper. 
Moreover, in terms of the long-time stability of such solutions it is clear that a low frequency analysis 
of the linearized operator is vital to obtaining suitable bounds on the corresponding solution operator of 
the nonlinear equation. 

To begin, notice that by Lemma [2] the linearized operator A has a T-periodic eigenvalue at the origin 
in the spectral plane of multiplicity (generically) three. Thus, as we allow small variations in the Floquet 
parameter, we expect that there will be three branches of continuous spectrum which bifurcate from the 
origin. According to Proposition [2 one of these branches must be confined to the imaginary axis, and 
hence will not contribute to any spectral instability. In order to determine if the other two branches 
bifurcate off the imaginary axis or not, we derive an asymptotic expansion of the function D{fj,,e^'^) 
for + |k| ^1. As a result, we will see that, to leading order, the local structure of the spectrum 
near the origin is governed by a homogeneous polynomial of degree three in the variables /i and n. We 
then evaluate the coefficients of the resulting polynomial in terms of Jacobians of various maps from the 
parameters (a, E, c) to the quantities T, A/, and P. To this end, we begin with the following Lemma. 

Lemma 4. // {T, Af, P}a,E,c 7^ 0, the equation D{fi, e'") = has the following normal form in a 
neighborhood of {fi, k) — (0, 0); 

_ (,«)3 + + tr ( M,,^(0)) - 5) - Af, P}^,E,c + + = (18) 

whose Newton diagram is depicted in Figure\^ 

Proof. Define functions a and 6 on a neighborhood of /i = by 

D{p, e*«) = -7^3 + {ai^l) - ■i)if + b{ti)r^ + D{^i, 1). (19) 
where rj = e*'^ — 1 is small. In particular, notice that 

a{fi) = tr(M(M)) 

bi^^) = i(tr((M(A.)-/)^)-tr(M(^)-/)2). 

^However, one should be aware of the recent works and [5] in which perturbation methods and more delicate 
functional analysis were used to analyze this problem for the gKdV in a way which provides more precise results than the 
above variational methods. 
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Figure 2: The Newton diagram corresponding to the asymptotic expansion of D{ii, e*'') = in a neighbor- 
hood of (/i, k) = (0,0) is shown to 0{\fj,\^ + k'^). Terms associated to open circles with dashed boundary 
are shown to vanish due to the natural symmetries inherent in ([7]) . The open circles with dark boundary 
are non- vanishing terms which are a part of the lower convex hull, and hence contribute to the dominant 
balance. The closed dark circles lie above the lower convex hull and thus do not contribute to the leading 
order asymptotics. See [3J and ^Sj for more details. 




Now, using the fact that the spectral problem ([7]) is invariant under the transformation {x, ^i) i— > (— x, — /i), 
it follows that the matrices M(— /i) and M(/i)^^ are similar for all /i e C. In particular, it follows that 

e"'"^ det(M(/^) - A) = -A^ det - i 

+ H-^^) - 3) - i) + K-m) (1 - i) + ^(-^' 

(A - l)-" - (a(-M) - 3) A (A - 1)^ + 6(-/i)A2 (A - 1) - X^D{-fi, 1) 

By comparing the O(A^) and ©(A^) terms above to those in (fT^ . we have the relations 

e-t'^/^aifi) = 2a(-/i) - b{-fi) - 3, 
_g-pT/c ^ -a{-fi) + b{-fi) - D{-^l, 1) + 2. 

Differentiating with respect to /i and evaluating at /i = immediately implies 6'(0) = and a'(0) = ^. 

Similarly, it follows that b"{Q) = a"(0) — — tr ( Mp^(O)) — The proof is now complete by Lemma 
Hand the fact that ?7 = + □ 

It follows that the structure of spec {A) in a neighborhood of the origin is, to leading order, determined 
by the above homogeneous polynomial in k and /i. Due to the triple root of _D(-, 1) at = the implicit 
function theorem fails, but can be trivially corrected by considering the appropriate change of variables. 
This leads us to the following theorem giving a modulational stability index for traveling wave solutions 
of d). 
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Theorem 3. With the above notation, define 





A = - ^tr(M^^(0))- (^-j j ^2tr(M^^(0))- (^-j ^ 

-27{r,Af,P}i^,, + 6{T,M,PK^„ ^3tr(M,,(0)) 5 

and suppose that {T, M, P}a^E,c 7^ 0. // A > 0, i/ien i/ie spectrum of the linearized operator A in a 
neighborhood of the origin consists of the imaginary axis with a triple covering. // A < 0, then cr(A) in 
a neighborhood of the origin consists of the imaginary axis with multiplicity one together with two curves 
which are tangent to lines through the origin. 

Proof. Since, to leading order in fi and n, the Evans function is homogeneous by Lemma 2] it seems 
natural to work with the projective coordinate y = Making such a change of variables, Lemma H] 
implies the equation D{fi, e'") = can be written as 




tr(M^^(0))~ - ] ~y^{T,M,P}a,E,c + i^E{K,y)^0 (20) 



, yT y^ 

where E[k,, y) is continuous in a neighborhood of the origin. Let 2/1,2,3 denote the three roots of the above 
cubic in y corresponding to E[k, y) = 0. Assuming A 7^ it follows that 2/1,2,3 are distinct and hence the 
implicit function theorem applies giving three distinct solutions of (|20p in a neighborhood of each of the 
2/1,2.3- In terms of the original variable fi, this gives three solution branches 

/^l,2,3 = -i2/i,2,3K + 0{k^). 

If A > 0, then 2/1.2,3 G M, giving three branches of spectrum emerging from the origin tangent to the 
imaginary axis. From the Hamiltonian symmetry of ([7]), the spectrum is symmetric with respect to 
reflections across the imaginary axis and hence A > implies these thee branches of spectrum must in 
fact lie on the imaginary axis, proving the existence of an interval of spectrum of multiplicity three on 
the imaginary axis. In the case A < 0, it follows that one of the roots, 2/1 say, is real while the other two 
2/2,3 occur in a complex conjugate pairs, giving one branch along the imaginary axis and two branches 
emerging from the origin tangent to lines through the origin with angel arg(— 12/2,3): see FigurelH □ 



Remark 3. The modulational instability index A derived above is considerably more complicated than 
the one derived by Bronski and Johnson JEI for generalized Korteweg-de Vries equation 

Ut = Uxxx - CUx + if{u))^ . (21) 

When considering the spectral stability of periodic traveling wave solutions of (j2ip . it was shown that there 
exists a modulational instability index AgKdV such that AgKdV < implies modulational instability, and 
AgKdV > implies modulational stability. In this case, the dominant balance is somewhat simpler due 
to the fact that the trace of the operator Mp(0) vanishes, and hence the K^/i term in the corresponding 
Newton diagram vanishes. As a result, the modulational instability index took the form 

AgKdV = ^ (tr(M^^(0)))' - i Q{T,A/,P},,s,c) 

where here D represents the Evans function for equation ([2T|l . It follows that the sign of D^^^{{), 1) does 
not effect the modulational instability of such periodic solutions of equation (I2ip . However, in the case of 
the gBBM equation, the fact that a'(0) 7^ seems to suggest that the sign of the sign of this orientation 
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Figure 3: When A < 0, the local normal form of spec(A) consists of a segment of the imaginary axis union 
with two straight lines making equal angles with the imaginary axis. Notice that these lines intersect at 
the origin, corresponding to the fact that D{fi, 1) = Od/xp). In this picture, the horizontal dashed line 
represents the real axis, and the dashed lines which are tangent to the straight dark lines represent the 
true spectrum of the linearized operator. 



index (see next section) has an impact on the modulational stability of the periodic traveling wave. In 
fact, in section 4 we prove exactly this fact in the case of a power-nonlinearity f(u) = u^~^^/{p+ 1): we 
prove the long-wavelegth periodic traveling wave solutions of ([1]) are modulationally unstable if and only 
if the limiting homoclinic orbit (solitary wave) is exponentially unstable. 

Our next goal is to use the integrable nature of ^ to express tr( M^^(O)) in terms of the underlying 
periodic traveling wave u. This is the content of the following lemma. 

Lemma 5. We have the following identity: 

i tr ( M^^(O)) = {T, P}e.,c + 2{M, P}a,E - V'{u.){T, M}a,E. 

Proof. This proof is essentially an extension of that of Lemma [2l Using the same notation, a straightfor- 
ward yet tedious calculation yields 

tr(M^^(0)) = tr (cof (<5 W(/i) W(0, 0)-')) 

= 2{T, P}e,c + 2{M, P}a,E + 2{r, M}a,c - -y'(U-){T, M}a,E 

c 

The proof is completed by noting that {T, M}a^c = {M, P}E,a since Tc = Pe and Mc = Pa- □ 

Thus, the modulational stability of a given periodic traveling wave solution to ([T]) can be determined 
from information about the underlying solution itself. Interestingly, from Lemma [5] the modulational 
instability index seems to depend directly on the turning point u_, a feature not seen in the corresponding 
analysis of the generalized KdV equation in [8J. However, one can trace this dependence back to the 
definition of wl in Lemma [5] and hence seems to be unavoidable. In the next section, we will analyze the 
above stability indices in neighborhoods of homoclinic orbits in phase space, complimenting the results 
of Gardner |12j by providing stability results in this limit. 
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6 Analysis of Stability Indices in the Solitary Wave Limit 

The goal of this section is to study the long wavelength asymptotics of the stability indices derived in the 
previous sectios. Throughout, we restrict ourselves to the case of a power-nonlinearity f{u) — u^^^ / (p+1): 
this restriction is vital to our calculation since in this case we gain an additional scaling symmetry. In 
particular, if v{x;a,E) satisfies the differential equation 

^---^+ (,+ l)W2) -"^-" + ^' 
then a straight forward calculation shows that we can express the periodic solution u(x] a, E, c) as 

1/2 „ p \ 



(23) 



This additional scaling allows explicit calculations of Pc, which ends up determining the stability of 
periodic traveling wave solutions of IlJ of sufficiently long wavelength. 

A reasonable guess would be that long-wavelength periodic traveling wave solutions of llj have the 
same stability properties as the limiting homoclinic orbit (solitary wave). However, as noted in the 
introduction, this is a highly singular limit and so it is not immediately clear whether such results are 
true. It is well known that the solitary wave is spectrally unstable if and only if p > 4 and 1 < c < co(p) 
for some critical wave speed co{p). It follows from the work of Gardner [12j that periodic waves below 
the homoclinic orbit in phase space which are sufficiently close to the homoclinic orbit are unstable if 
the solitary wave is unstable. In particular, it is proved that the linearized operator A for the periodic 
traveling wave u with sufficiently long wavelength has a "loop" of spectrum in the neighborhood of any 
unstable eigenvalues of the limiting solitary wave. Thus, Gardner's analysis deduces instability of long 
wavelength periodic waves from instability of the limiting solitary waves. The results of this section 
compliment this theory by also proving that the stability of the limiting wave is inherited by nearby 
periodic waves below the homoclinic orbit in phase space. 

In terms of the finite- wavelength instability index, it seems reasonable by Theorem[T]to expect that for 
periodic traveling wave solutions below the seperatrix of sufficiently long wavelength, {T, Af, P}a,E,c < 
for if and only if p > 4 and 1 < c < co{p). What is unclear is whether such a result should be true 
for the modulational instability index A. Indeed, although Gardner's results prove that the spectrum of 
the linearization about a periodic traveling wave of sufficiently long wavelength in the neighborhood of 
the origin contains the image of a continuous map of the unit circle, to our knowledge it has never been 
proved that this map is injective. Thus, it is not clear from Gardner's results whether a modulational 
instability will arise from this eigenvalue since it is possible this "loop" is confined to the imaginary axis: 
we show that in fact one has modulational instability in this limit precisely when the limiting solitary 
wave is unstable. 

The main result for this section is the following theorem, which is based on asymptotic estimates of 
the instability indices derived in section 3. In particular, we prove the sign of both instability indices 
in the solitary wave limit is determined by the sign of -^^P, where P — P{a, E, c) is the momentum 
of the periodic wave u{x;a,E,c). The proof is based on a more technical lemma, which shows that 
-^M{a, E,c) < for waves of sufficiently long wavelength, i.e. for a,E sufficiently close to zero. We 
begin by outlining the proof of the following theorem, and then fill in the necessary lemma's afterward. 

Theorem 4. Let f{u) — uP^^ for some p > 1 and let u{x] a, E, c), with a sufficiently small, be a periodic 
solution of the traveling wave ODE ^ which corresponds to an orbit below the homoclinic orbit in phase 
space. IfT{a, E, c) is sufficiently large, then u is modulationally and nonlinearly stable for all 1 < p < 4. 
Forp > 4, there exists a critical wave speed c{p) > 1 such that ifT(a, E, c) is sufficiently large, the solution 
u is orbitally and modulationally stable if O c{p), while it is modulationally unstable for 1 < c < c{jp). 
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Proof. Without loss of generality, suppose that the periodic solution u{x; a, E, c) of ^ corresponds to a 
branch cut of the function \/ E — V{u\ a, c) with positive right end point. Now, when a and E are small 
there are two turning points r\ , r-i in the neighborhood of the origin and a third turning point ra which 
is bounded away from the origin. In the solitary wave limit a, £^ ^ a straight forward calculation gives 



that r2 — ri = \ \J <^ — 2(c — From this, it follows that the period satisfies the asymptotic relation 

T(a,E,c)=0{\^{a^ -2{c-\)E)) . (24) 

To see this, notice that by our assumptions on ri, r2, and r-^ we can write E — Y{u\ a, 1) = (r — ri)(r — 
''2)('''3 — ^^Q{t) where Q(r) is positive on the set [ri,r3]. The period can then be expressed as 



T{a,E,c) = \f2c \ 

J To 



dr 



^J{r - ri)(r - r2)(r3 - r)Q(r) 



T2-ri \Jr{r - (r2 - ri)){rz - ri - r)Q{r + ri) 

'■^-'^1 dr 

i+<5 \/r{r - (r2 - ri))(r3 - ri - r)Q{r + ri) 



The integral over the set (r2 — ri + (5, — ri) is clearly 0(1) as (a, E) — * (0, 0). For the other integral, 
notice that 

(ra - ri - r)Q{r + ri) > 
on the set [r2 — ri,r2 ^ ri + S]. Thus, in the limit as (a, E) oo we have 



T{a,E,c)^ / 



^yr{r - (r2 - ri)) 



= -41n(4(r2 - ri)) + 2 ln(2(v/r2 - ri + 5 + %/^)) 
~ -ln(r2 - ri) 

from which (jM]) follows. Similar computations yield the following asymptotic relations for a and i? 
sufficiently small: 

F(a,i?,c) = 0(l) 

M{a,E,c) = O (aln (a^ - 2(c- 1)^;)) 
Ma{a,E,c)^0 





a2 


- 2(c- 


l)i? 




a 




a2 


- 2(c- 


l)i? 




1 




a2 


- 2(c- 


1)E 




E 




a2 


- 2(c- 


1)E 



Ta{a, E,c)^0[ ) = Msia, E, c) 

TE{a,E,c) = 



T,{a,E,c) = O ( ]=2Pe + -T. 

Thanks to the above scaling we know Mc = 2Pa + ir can be expressed as a linear combination of M, 
Ma, and Me — Ta- Similarly, Pc can be expressed as a linear combination of P, 2Pa — Mc ^ \T and Pe- 
It follows that the asymptotically largest minor of {T, Af, P}a,E,c for a and E small is —TEMaPc and, 
moreover, 

tr(M,,^(0))-rBP,. 
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In Lemma [S] below, we will show that Ma < for such a and E. Moreover, it is clear that for such 
periodic waves with sufficiently long wavelength satisfy Te > 0. Therefore, it follows that both stability 
indices are determined by the sign of Pc(a, E, c) in the solitary wave limit. The theorem now follows by 
Lemma [7] below. □ 

Remark 4. Notice that the fact that the finite wavelength instability index is determined by the sign of 
Pc in the solitary wave limit is not surprising, since this is exactly what detects the stability of the limiting 
solitary waves. What is surprising is that the same quantity controls the modulational stability index in 
the same limit. As mentioned above, it has not been know if the instability of the limiting solitary wave 
forces a modulational instability: the answer is shown to be affirmative by Theorem^ 

In order to complete the proof Theorem |31 we must prove a few more technical lemmas. The first is 
used in showing that the sign of the modulational and finite- wavelength instability indices are determined 
completely by the sign of Pc{a, E, c) in the limit as a, E tend to zero. This is the result of the following 
lemma. 

Lemma 6. For a, E sufficiently small, Ma{a, E,c) < for all c > 1. 

Proof. Notice it is sufficient to prove ^M(a,0, c) < for all c > 1 and a sufficiently small. Now, 
M(a, 0, c) can be written as 



M(a, 0, c) = V2c 



^^"'■'^^ ^/u du 



2 " (p+l)(p+2)" 



where r(a, c) is the smallest positive root of the polynomial equation a + ^2^r — (p-|-i)'(p+2) t^^^ ~ 0. 
Setting a= {{p+ l){p + 2))i/Pa, we have 

fr(a,c) r- , 

M {a, E, c) = ((p + l){p + 2)f/PV2'c 



c-1 

^ 2 

speed c > 1, r{a,c) is a smooth function of a for a sufficiently small and satisfies 



where r{a, c) is the smallest positive root of the polynomial a + ^^r — rP'^^ — 0. Notice for a fixed wave 



1\'/P 2a 



The goal is to now rewrite the above integral over a fixed domain and show that the integrand is a 
decreasing function of a for a fixed c > 1. 

Making the substitution u —^ r{a, c)u yields the expression 

M(a,0, c) _ ^/udu 



{{p+l){p + 2)Y/'P^c Jo ?■(«, c)-3 + (^) u ¥{a, c)-^ - r{a, c)p-^ uP+^ 
Now, we can use the above expansion of r(a, c) to conclude that 

^ (^a f(a, c)-3 + u ria, c)-' - r{a, c)p-' uP+'^ 

is positive on the open interval (0, 1) for all p > 1 and c > 1, which completes the proof. □ 

Finally, it is left to analyze the asymptotic behavior for a fixed wavespeed c > 1 of the quantity 
Pc(a, E, c) in the solitary wave limit. This is the content of the following lemma. 
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Lemma 7. In the case of power non-linearity f{u) = u^^^ /{p + 1), we the momentum P — P{a,E,c) 
satisfies 

1)2/^-1/2,1/2^ Uc + p){c-l)p\ 

-P{a,E,c) = — [4c-p+ ^ +0{\a\ + \E\) 



dc ' ' 2pc{c - 1) V (4 + P)c 

in the solitary wave limit {a,E) (0,0), where I{r) — J^^sech^{x)dx. In particular, for a and E 
sufficiently small, if p < 4 then -^P {a,E,c) > for all c > 1 while if p > 4 then -^P {a,E,c) < for 
1 < c < cq{p) and -^P (u; a, E,c) > for c > Cf){p), where 



Cob) = 



4 + 2p 



Proof. The proof is based on scaling and a limiting argument, as well as a modification of the analysis 
in [IHj. To begin, let v = v{x;a,E) satisfy the differential equation ([2H) so that u{x;a,E,c) can be 
expressed via scaling as in (j23p . and assume with out loss of generality that a; = be an absolute max of 
v{x; a, E, c). Clearly, the solitary wave limit corresponds to taking (a, E) — > (0, 0) with fixed wave speed 
c > 1. Notice that on any compact subset F of K, we have 

vix; a, E) ^ {^l±3k±}l^ ^ech^/. + 

uniformly as {a,E) (0,0) on F for some xq E M. Using ([22|) . it follows that 

r.T((c~l)/c)i/^ 



/ v'^{x)dx = V2 / 

Jo Ju- 



where u± are the roots of £' + — (p_|_i)'(p_|_2) v^^^ + av — satisfying the original hypothesis of the 
roots u± of E — V{u;a,c) = 0. Since /5t('(c-iVc)i/V2 ^^^"^^^^ ~ ^(1) («;^) ~^ (0,0), the dominated 

convergence theorem along with the fact that j_T({c-i)/cY/y2 '^'^{^)dx is a function of a and E implies 
that 

T((c-i)/c)i/V2 f{n + 2){v+l)V^ I /4\ 

^T((c-l)/c)i/V2 V 2 / P VP/ 



(p + 2)(p+l) y/^ 1 /4 



z;,(a;; a, i?)^da; = '-i^ _— / - + 0{\a\ + |£;|) 



Similarly, it follows that 

rT((c-l)/c)i/V2 

^T((c-i)/c)i/2/2 \^ 2 y 4+p 

Using we now have 

(p + )(p+ ) \ /" (^u{x\a,E,c)'^ + U3;(x;a,E,c)'^^dx 

T/2 



: (C - l)2/p-l/2cl/2j f 1^ + (c ^ i)2/p+l/2,-l/2^ ^ f 

0(|a| + |i?|) 



Stability of Periodic gBBM Waves 



23 



&s a,E ^ 0, and hence it follows by differentiation that 



-T/2 

(c-1 



\2/p-l/2„l/2 7- ^ 

^ VpJ f^^_p , (4c + p)(c-l)p 



2pc(c-l) V (4 + p)c 

+ 0(|a| + |£;|) 

as claimed, where we have used that TcU^ — Odal + I-Bl). The lemma now follows by solving the quadratic 
equation (4 + p)c(4c — p) + (4c + p)(c — l)p = for c and recalling the restriction that c > 1. □ 

The proof of Theorem|3]is now complete by Lemmas El and [71 As a consequence, the finite-wavelength 
instability index {T, M, P}a.E.c seems to be a somewhat natural generalization of the solitary wave 
stability index, in the sense that the well known stability properties of solitary waves are recovered 
in a long-wavelength limit. Moreover, this gives an extension of the results of Gardner in the case of 
generalized BBM equation by proving the marginally stable eigenvalue of the solitary wave at the origin 
contributes to modulational instabilities of nearby periodic waves when ever the solitary wave is unstable. 

Finally, we wish to make an interesting comment concerning the elliptic function solutions of the 
modified BBM (mBBM) equation with f{u) = u^. In the special case in our theory of a = 0, it turns out 
that all periodic traveling wave solutions u{x; 0, E, c) can be expressed in terms of Jacobi elliptic functions: 
when > such solutions correspond to cnoidal waves, expressible in terms of the cn(x, k) function, 
while when E < such solutions correspond to dnoidal waves, expressible in terms of the function 
dn{x,k). By Theorem HI it follows that the dnoidal waves with 1 — k sufficiently small(corresponding 
toO<— ii^<Cl) are always modulationally stable and nonlinear ly (orbitally) stable to perturbations 
with the same periodic structure. What is unclear from our above analysis is whether the cnoidal waves 
with k sufficiently close to 1 (corresponding to < i? ^ 1) are exhibit a sense of stability: such waves 
correspond to periodic orbits outside the seperatrix in phase space, and hence Theorem 2] does not apply 
to such solutions. Moreover, such waves do not approach any particular solitary wave in any uniform 
way, even on compact subsets, and hence it is not clear if the stability of the limiting solitary waves at all 
influences the cnoidal waves. However, notice that for k sufficiently close to unity one must have Te < 
and hence by the analysis in the proof of Theorem [4] we have A < 0. Therefore, the cnoidal waves of 
the mBBM with sufficiently long wavelength are modulationally unstable. Similarly, in this solitary wave 
limit one has 

lim sign({r,M,P},,£,,(0,£;,c)) = lim sign (A4(0, c)) 

and hence one has spectral instability to periodic perturbations if Ma{0, E,c) < for E > sufficiently 
small. Using the complex analytic calculations of Bronski, Johnson, and Kapitula [5], we can solve the 
corresponding Picard-Fuchs system and see that 

disc{R{u;0,E,c))Ma{0,E,c) _ ^ 

(c- l)T(0,£:,c) ■ ^ ' 

See the appendix for more details of this calculation. Since disc(i?(u; 0, E, c)) < for such solutions, we 
find that Ma < 0. Thus, the cnoidal wave solutions of the mBBM equation of sufficiently large wavelength 
are exponentially unstable to perturbations with the same periodic structure. In particular, it seems like 
the index sign(T£;) serves as Maslov index in this problem: when Te > it is Pc < which signals the 
stability near the solitary wave, while it is Pc > wheiH Te < 0. This index arises naturally in the 



^While we have not actually proven this here (we only showed this is the case for cnoidal solutions of mBBM), we suspect 
this is the case for all nonlinearities of the form /(«) = n^'°+^ for fc 6 N. For examples where this phenomenon is studied 
in more cases see [9]. 
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solitary wave setting: see |10| for details and discussion. In particular, notice that in the solitary wave 
setting of Pego and Weinstein and in the long wavelength analysis of Gardner, the Maslov index always 
has the same sign: however, as seen here this is not the case when considering the full family of periodic 
traveling waves of the gBBM. 

7 Concluding Remarks 

In this paper, we considered the stability of periodic traveling wave solutions of the generalized Benjamin- 
Bona-Mahony equation with respect to periodic and localized perturbations. Two stability indices were 
introduced of the full four parameter family of periodic traveling waves. The first, which is given by 
the Jacobian of the map between the constants of integration of the traveling wave ordinary differential 
equation and the conserved quantities of the partial differential equation restricted to the class of periodic 
traveling wave solutions, serves to count (modulo 2) the number of periodic eigenvalues along the real axis. 
This is, in some sense, a natural generalization of the analogous calculation for the solitary wave solutions, 
and reduces to this is the solitary wave limit. The second index, which arises as the discriminant of a 
cubic which governs the normal form of the linearized operator in a neighborhood of the origin, can also 
be expressed in terms of the conserved quantities of the partial differential equation and their derivatives 
with respect to the constants of integration of the ordinary differential equation. This discriminant detects 
modulational instabilities of the underlying periodic wave, i.e. bands of spectrum off of the imaginary 
axis in an arbitrary neighborhood of the origin. 

In our calculations, we heavily used the fact that the ordinary differential equation defining the 
traveling wave solutions of ([T|) have sufficient first integrals, and as such is doubtlessly related to the 
multi-symplectic formalism of Bridges [7]. Many of the ideas from this paper were recently used by 
the author, in collaboration with Jared C. Bronski, in the analogous study of the spectral stability of 
periodic traveling wave solutions of the generalized Korteweg-de Vries equation [5]. Thus, the general 
technique of this paper has been successful in two different cases. These successes gives an indication 
that a rather general modulational stability theory can be developed using these techniques. It would 
be quite interesting to compare the rigorous results from this paper to the formal modulational stability 
predictions of Whitham's modulation theory [5T]. Moreover, the techniques in this paper could eventually 
lead to a rigorous justification of the Whitham modulation equations for such dispersive equations. 

It is not known if our techniques can be extended to equations which are deficient in the number of 
first integrals. An example would be the generalized regularized Boussinesq equation 

Utt - Uxx - {I{u))^^ - Uxxtt = 0. 

A straight forward calculation shows that although this equation admits a five parameter family of 
traveling wave solutions with a Jour parameter submanifold of periodic solutions. Although one can 
still build a four dimensional basis of null-space of the corresponding linearized operator, it is not clear 
whether the resulting stability indices can be directly related back to geometric information concerning 
the original periodic traveling wave. Such a result would be very interesting, as it would allow the use of 
such techniques to several more classes of partial differential equations admitting traveling wave solutions. 

Acknowledgements: The author would like to thank Jared C. Bronski for his many useful conver- 
sations in the early stages of this work. Also, the author gratefully acknowledges support from a National 
Science Foundation Postdoctoral Fellowship under grant DMS-0902192. 

8 Appendix 

In this appendix, we outline the general complex analytic methods used to derive the identity in equation 
()25p . To begin, consider ([T]) with f{u) ~ u"^ and let u{x]a,E,c) be a given periodic traveling wave. To 
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begin, we notice that that in this case T, M, and P can be represented integrals on a Riemann surface 
corresponding to the classicaUy allowed region for \J E — a, c), where V(u; a, c) = i? + au+ (^^) — 
In particular, defining R{u; a, E, c) := E — V{u; a, c) we have the following integral representations: 

T{a,E,c) 
M{a,E,c) 




P{a,E,c)^\ 

where K is the classical action defined in ([6]) and the integrals are taken around the classically allowed 
region corresponding to the periodic orbit a, i?, c). By differentiating the above relations with respect 
to the traveling wave parameters, we have the following identities: 

/ 1 fc [ u du ^ fc I du 1 fc [ u^ du 1 

„ , ^ ^ / 1 fc f du ^ 1 fc f du 1 

In particular, we see that each of the above derivatives can be expressed as linear combinations of the 
classical action K, the /c*'* moment of the solution u 



fJ-k 

and an integral of the form 

h = 



/u^ du 



u^ du 



R{u)y^' 

We now show that in the case of the mBBM equation, we can find a closed system of seven equations for 
the integrals Ik in terms of the traveling wave parameters and the quantities T, M, P and K. To this 
end, first notice that for any integer j > I one has the identity 

/u^ du f u^ R{u)du 

= E Ij+a Ij+i + {—^^ h+2 - 

Moreover, integration by parts yields the relation 

/ u'-^ du f uiR'(u)du 

= alj + (c- l)/j+i -Ij+2- 
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This yields a linear system of seven equations in seven unknowns {Ij}j^o known as the Picard-Fuchs 
system: 

2M 



/ E 
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c-l 
2 
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4 





^ \ 




( ^° 
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The matrix which arises in the above linear system is known as the Sylvester matrix of the polynomials 
R{u) and R'{u). In commutative algebra, a standard result is that the determinant of the Sylvester matrix 
of two polynomials P{u) and Q{u) vanishes if and only if they have a common root. If we now restrict 
ourselves to a = 0, it follows that the above matrix is invertible for all periodic traveling wave solutions 
u{x; 0, E, c) of the corresponding mBBM equation. In particular, we can solve the above Picard-Fuchs 
system and use the fact that Ma = — \-\/%l2 to see that 

disc:(/?((/:().£,c))il/,,(0,S,c) , 
lim — '■ — _ — — r- = 1 



(c-i)r(o,£;,c) 



as claimed. 
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